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ABSTRACT 


The problem of applying the theory of systems analysis 
to the marine propulsion plant is investigated in this thesis. 
It is believed by the author that the use of this theory can 
be of considerable assistance to the naval engineer in giving 
him an understanding of the fundamental behavior of the system 
and will provide a scientific basis for design criteria and 
operational policies. 


There are five basic areas of investigation covered in 
the thesis which all contribute to the organization of the 
overall method of accomplishing an operational analysis of a 
complex mechanical system. The following is a concise digest 
of each of the sections. 


1. The derivations of the basic mathematical formulas 
for reliability determination are given. The concept of com- 
ponent failure distribution is presented and the exponential, 
normal, and Weibull failure models are discussed. 


2, The importance of defining the objectives of the 
analysis clearly and completely is shown, The interpretation 
of component operating conditions and the resolving of 
suecess and failure standards are explained and an outline 
for the formulation of the necessary input data required to 
undertake the mathematical part of the analysis is offered. 

A hypothetical combined diesel-gas turbine propulsion plant 
is used to illustrate the proper procedures for performing a 
reliability analysis. 


3, The mathematical procedures for determining 
system reliability by the correlation of the individual 
system components behavior are explained. In addition to 
presenting methods for solving serial, parallel, and Bayesian 
systems, formulas are derived for the reliability analysis of 
periodically operating components and of standby units. 
Representative values are postulated for the diesel-gas 
turbine plant and reliability curves are plotted which demon- 
strate the effect of component wear-out, reliability as a 
function of time, and the individual effect of specific com- 
ponents on the overall system reliability. 
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4, Two techniques for ascertaining maintenance 
schedules are formulated. The first technique yields the 
optimum time to perform preventive maintenance in order to 
maximize system availability or minimize expected repair 
costs. The second procedure deduces by dynamic programming 
an optimal decision policy on whether or not to repair the 
system dependent upon the number of trips made since the last 
repair and the number of trips still to be accomplished. 


5. A brief investigation is made into the problem 
of components subject to the interaction of their failure 
distributions. Four possible interaction functions of 
instant failure rate are proposed and analyzed as to their 
effect on reliability and mean time to failure. An example 
is used to show how important this interaction effect can be 
when it is time dependent,such as would be caused by gradual 
wear. 


In conclusion this thesis presents a procedure for ana- 
lyzing a complex mechanical system using hypothetical 
examples to better illustrate the utility of the results to 
the naval engineer. | 


Thesis Supervisor: Ernst G. Frankel 
Title: Assistant Professor of Naval Architecture 
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NOMENCLATURE 


Meaning 
Reliability 


Instant failure rate 

Probability density function of 
failure 

Probability distribution of 
failure 

Mean time before failure 

Time 

Exponential instant failure rate 


Constant 


specific nomenclature for diesel-gas turbine reliability analysis 


p(t) for diesel engine 

p(t) for gas turbine 

о(%) for heat exchange 

p(t) for reduction gear 

p(t) for ibe o11 pump 

p(t) for standby lube oil pump 

p(t) for fresh water pump 

p(t) for salt water pump 

Mean time of non use of the gas 
turbine 

Probability that the gas turbine 
will not be in use at time t 
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INTRODUCTION 


Operational Analysis and Reliability Engineering 
им During and shortly after the Second World War, due to 


the growing complexity of all types of operations in 
general, a branch of mathematics termed operations analysis 
became prominent. Probability theos? forms the basis for 
this field of science. If a definition is desired, it may 
be said the operations analysis is a scientific study of 
operations to better understand their behavior por the 
twofold purpose of predicting nei future results due to 
changes in the system and controlling the operation to im- 
prove its result. This thesis is concerned with the use of 
operations analysis methods to solve problems in the field 
of reliability engineering, 

| In any system, mechanical — one of the most 
important parameters determining its overall value is the 
reliability that the system can be expected to have during 
its operation. Although the exact definition of reliability 
is usually dependent upon who is the definer, the accepted 
definition endorsed by the American Soctety for Quality 
Control [1] states that reliability is the probability of a 
device operating within specified limits for the time and 
operating conditions imposed on the device. Operations 
analysis techniques are used in reliability engineering for 
the purpose of predicting reliability and optimizing the 
system from a reliability standpoint. While reliability 


analysis is a relatively new subject which has generally 
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only been applied to the electronics and missile field, it 
has the potential of being extremely useful in the under- 
standing and improving of all types of mechanical systems 


including marine propulsion plants. 


Importance of Reliability Analysis 
The ability to be able to predict reliability is tanta- 


mount to understanding the underlying causes of unrelia- 
bility. The ИРЕТИ of unreliability are numerous; the 
most obvious ones being cost, lost time, and danger to 
personnel. The cost involved not only is that of the device 
to be replaced or repaired, but also the cost of the mainte- 
nance in material, manpower, and training with the added 
expense of having to continuously monitor the equipment for 
possible failures. The value of the analysis is that in 
choosing some specific system to perform a task, a relia- 
bility analysis of the system will furnish data on which to 
base a decision. For instance, the Armed Services mainly 
places the emphasis for selecting a component or system on 
the low bid price. Normally, reliability is only given 
token consideration in the decision. The question srises as 
to whether a reliability investigation of the device should 
be made since isn't it possible that it would be aqya pie 
to pay more initially for a reliable system that in the long 
run will give better availability and less maintenance costs 
than a cheap unreliable system? Conversely, the possibility 
arises as to whether in some applications too much relia- 


bility is required of a device. Certainly the adage that 
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"a chain is only as strong as its weakest link" applies in 
a system and, for example, designing a valve with a safety 
factor of 5 when the piping only has a safety factor of 2 
is a waste of time and money. It is therefore evident that 
reliability apportionment should enter into any preliminary 
system design. 

Another problem inherent to this day and age is that of 
complexity. The major problems to be solved have become 
more intricate as time passes „ра it is заг natural that 
there is an increase in the complexity of the equipment 
designed to meet these problems. The more complex equipment 
is, the greater number of components are required, and the 
more likely it is that one of these components will fail, 
Thus, an increase in complexity leads to a decrease in 
reliability. Of what value is complex equipment with theo- 
retically high performance if failures of its components 
keep the piece under repair most of the біте, It can 
therefore be seen that engineering improvements for relia- 
bility will in many cases improve operability, maintaina- 
bility and productivity with an overall net gain in system 


worth. 


Other System Analysis Results 
Results other than reliability may also be obtained from 


an operations analysis of a system. The question of main- 
tainability of the system showld normally be investigated 
unless for some special reason the system cannot be main- 


tained during its operating cycle, i.e. missiles, torpedoes, 
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etc. A study of maintainability will give answers to such 
questions as; 

l. Is preventive maintenance worthwhile for the 
particular system? 

2, Lf worthwhile, at what time intervals should it be 
scheduled? 

3. What is the optimum efficiency that can be expected 
from a system preventively maintained? 

in conjunction with preventive maintenance scheduling, 
the renewal theory of operations analysis [2] will give the 
expected number of replacements for a component over its 
lifetime, t. This gives a good basis for making an estimate 
of the number of spare parts to be stocked. 

Some other results that are obtainable are; 

1. The effect of redundancy on reliability. 

о. The system availability that can be expected. 

3, Optimization procedures for reasons of reliability 
or cost. 

4, Confidence levels that can be applied to the 
results of the analysis. 

5, Information for the determination of equipment 
quality control specifications. 

It is seen that by the methods of operations analysis, 
a deep insight into the behavior of a system can be obtained 
which will produce not only a better understanding of the 
system but also give a firm basis for making decisions to 


improve its overall value. 
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Present State of the Art 

At the present time the state of the art of operational 
analysis of mechanical systems is very sketchy. Many mathe- 
matical models have been devised of hypothetical systems 
which have been analyzed thoroughly but the analysis of real 
systems has only been accomplished in the electronics and 
missile fields and even there only in limited situations. 
There are two main reasons for this lack of progress. First 
and foremost is the fact that the value of systems analysis 
and especially reliability engineering has just in the past 
few years come into prominence. Now that the mathematical 
tools have been developed and the value realized, the majori- 
ty of large companies are beginning to organize reliability 
analysis departments, The second factor and possibly the 
hardest to solve is the cost and time involved in running 
the testing programs necessary for the analysis. In order 
to obtain any results from an analysis, failure tests must 
be run for the devices under investigation and reams of data 
must be evaluated, Although this statistical aspect of 
operations Вос 5516 will not be discussed in this thesis, 
it should be mentioned that there are sampling procedures 
that can be used effectively to make data correlation ef- 
ficient, Inherently, the cost of initiating a reliability 
program will be high but just as inherently is the fact that 
in the long run, the program will more than pay for itself 


in its results. 
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Thesis Objective 
The object of this thesis is to show how the methods of 


operations analysis can be employed in the investigation of 
a marine propulsion plant and of what worth such an investi- 
gation will be. The overall study of this question will be 
divided into five parts; 

1, Reliability analysis basic concepts. 

2. Formulation of the problem of determining system 
behavior by correlation of individual system component 
behavior. 

3. Mathematical determination of system reliability. 

4. Effect of preventive maintenance, 


5. The problem of interacting components. 
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PROCEDURE 


I. Reliability Analysis Basic Concepts 


Before introducing the method by which to solve the 
reliability analysis of a complex propulsion system, a 
short explanation of the theory behind the mathematical 
derivation of reliability will be given. This is deemed 
necessary since reliability analysis is essentially a new 
field and although it has as a foundation the mathematics 
of probability and of statistics, there are some intrinsic 
differences. The standard reliability terms are defined in 
the Nomenclature section of this thesis, however since the 
basic mathematical formulation is short, it will be included 
in the text. 

To determine reliability mathematically, the following 
is presented; 

Let 

x = P(Xy 2X52 0-05) be a vector of performance 
| | characteristics of the 
component 
А - region of the sample space representing the 
satisfactory performance of the component. 
w= £(W, Wos 00.» W,) be vectors specifying the 
| upper and lower limits оп 
2 f (Z4 s Zo» eo o3 Zy) satisfactory performance 

t = time of the interval under investigation 
f(x,t) = probability density function of the performance 

characteristics 


R = reliahilitv of the component 
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therefore 
R(x C A; t) = Prob.(w < х zZ) CUN 
= A 


employing the proper definitions [20] 


- 


R(t) = f" f(x) ax = 1 - H f(x) dx =1 - F(t) 
furthermore -fË p(z) dx 
R(t) = e ы | 
where | 
p(t) = i = - 5 (1a R(t). 
-[* p(x) ах 
S) -p(E)e 


The other parameter of reliability mathematics that is 
of great importance is the components! mean time before 
failure (MTBF). 
ти T. = mean time before failure 

Е = expected value 
then 
Ta = E(t) - ^" t f(t) dt = Г t dR(t) = Г R(t) dt 


Im the event that the component has been operating for a 
period previous to the time of investigation; 


Defining t = O time component is new, or system starts 


operating 
t = t time of investigation 
$ = time to point of interest 
R(t/t. ) = reliability of the component at time t given 


the component has operated to time бү. 


thus 
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R(t/t,) = REST 


2 R(t 
г t (6-61) Же г dt 


RET г R(t) dt 


1 


and 


= 
lI 


Component Failure Distributions 

From the preceding paragraphs it can be seen that relia- 
bility and mean time to failure can be evaluated from the 
probability density life distribution of the component, In 
Ева following paragraphs some of the more important mathe- 
matical approximations of component life distributions are 
discussed. These distributions have been derived from ex- 
tensive analysis of operating data and life tests. It must 
be remembered, however, that they are just statistical models 
characterizing a physical phenomenon and are, at least, a 
good approximation of the reality. For a better understanding 
of the real situation, confidence limits are associated with 
each of these approximations. 

Normally, the characteristic parameter used to describe 
life distributions is the instant failure rate, p(t). The 
time dependence of p(t) that is representative of a normal 


component is shown in Fig. I. 
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FIGURE 1 


Representative p(t) vs t plot 


o 


III - 





early failure region which results from manufacturing 
flaws, in-transit damage, etc, which cause normal 
failure at an early stage in the operating cycle. 
Usually in real situations there is a "debugging" or 
run-in period which eliminates this problem. 

chance failure region caused by random failures which 
are normally considered independent of time and 
having a fixed probability of occuring at any time in 
this region. 

wear-out region due to the gradual deterioration of 


the component with time. 


Exponential Distribution 


This is the most widely used life distribution in 
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reliability analyses. The reasons for this are it; 

1. characterizes components that have been optimized 
to limit failures, 

2, characterizes components having a predominence of 
human errors; 

3. characterizes complex components, 

|l. characterizes components consisting of parts of 
mixed ages, 

5. approximates what is usually the behavior in the 
chance failure region, 

6, is the easiest mathematically to work with. 


The equations associated with the exponential distribution 


are; 
p(t) = constant = c 
f(t) = p(t) erplt)t 
R(t) = eet 
E- 597 
FIGURE II 


Exponential Distribution 
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Normal Distribution 


This distribution is generally employed to describe the 


wear-out phenomenon. It also has been found to apply to; 
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1. components exhibiting homogeneous deterioration 
properties, 

2, components subjected to small variations in environ- 
mental severity, 

3, components whose failures occur at times well 
removed from t = O, and whose MTBF is large 
compared to its standard deviation. | 

The equations associated with the normal (Gaussian) distri- 


bution are; 








M = mean lifetime 
o = standard deviation of M 
20 
е 
р(%) ---------5--- 
- (т-м)2 
Ге = дт 
t 
В (t-M)* 
1 Po 
f(t) =— е 
2 аУ?т 
„ AZM С 
R(t) = s" Ге =o dt 
| о42т t 
Tp =M 
2 px as 
if erf(x) = Ld e du 
then 
1 $-М 
R(t) = = [1 - erf(——)] 
2 oN? 
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FIGURE III 
Normal Distribution 
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t t t 
Weibull Distribution : 

Although this distribution has rarely been used in 
reliability analyses to date, statistical investigation of 
life tests of mechanical tests shows that this distribution 
seems to be applicable in many cases. The distribution can 
be described by a two parameter method or a three parameter 
method. It includes as a special case the exponential 
distribution and can approximate a normal distribution 
(в - 3.5). | | 
For the three parameter case; 

where а - Scale parameter t > Y, G, P > O 


P 
y 


slope parameter 


location parameter 


Qa . 


0) = ВО?" 
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_ (s-y)Ë 


T 
R(t) =e 


т. = ү + а? г(8+1) 


FIGURE IV 
Weibull Distribution 
FW) ee КИ 
К 
el 
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Other Distributions 

Some other distributions occasionally used in relia- 
bility studies are; 

1. lognormal 

2, gamma or beta 

3. mixed Weibull 
however, the exponential, normal, and Weibull are the most 
important of the distributions and are the ones used for 
examples in this thesis. 
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11 Formulation of the Problem 
Defining the Problem 

One of the most critical parts of performing an analysis 
of any system is the complete definition of the problem to be 
studied. This definition basically must answer three 
questions; 

1. What are the objectives of the analysis? 

2. What components are to be considered comprising the 
system? 

3. What will be considered the failure events of the 
system and components? 

These must be clearly and exactly defined if clear and 
exact answers are to be expected. In some cases if just 
figure of merit or trend results are desired, greater liber- 
ties may be taken in the definitions, Once these three 
questions are answered, the problem has been defined and the 
results of the investigation are completely dependent upon | 
these definitions. 

Objectives 

When defining the objective of the analysis, the question 
is, what is wanted and what should be expected from the 
analysis. In order to better illustrate how to accomplish an 
operations analysis of a system, a hypothetical marine pro- 
pulsion plant will be postulated and investigated. It must 
be realized by the reader that this example will be necessarily 
simplified since the object is not numerical results but 
instead the method of the solution. 
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The objective of this analysis will be to find the 
reliability of a combined diesel and gas turbine marine 
propulsion plant, The result desired is how does the 
reliability of the overall plant vary with time in орег- 
ation and how should one go about improving this relia- 
bility. Naturally in this example the value of the 
numerical results. are useless since there is no statisti- 
cal basis for the numbers used. In an actual analysis, 
statistical life tests would have to be performed on all 
the components involved. From these statistical tests, 
confidence limits can be found and applied to the final 
results. Methods of finding confidence levels are dis- 
cussed in Lloyd and Lipow [2]. 

System | 

The next step in the problem is to define the system, 
The system can be defined in as basic or complex terms as 
the investigator deems necessary in order to obtain the 
results desired. In this example the system will be com- 
prised of the following components; 

і, diesel engine 

o lube oil pump 
ə fresh water pump 


. salt water pump 


ə gas turbine 


2 
3 
4 
5, heat exchanger 
6 
T. reduction gear 
8 


. Standby lube oil pump. 
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From the definition of reliability it is seen that in 
defining the system, the operating conditions which are 
within the prescribed limits must also be enumerated. 

There are five basic operating conditions that should be 
taken into accounts; 

1. environmental 

2, conditions imposed by the operator or user of the 

system 

3. functional dependence of the components within the 

defined system 

4. interaction of component failure distributions 

5. maintenance policies. 

It is altogether possible that more than one combi- 
nation of these conditions may want to be investigated. Іп 
this case different reliabilities can be expected for the 
different combinations. For the purpose of this investi- 
gation the following operating conditions are defined; 

Environmental: The only extraordinary environmental 

conditions imposed on the system will be the possi- 

bility of shock with the probability of Р of 
occurring. 

Imposed conditions: The system must be able to 

operate for 20% of its total operating cycle above 

cruising power. 

Functional dependence: The functional dependence of 

system operation upon the components is; 

1. Diesel must perform successfully for the system to 


operate. 


о Rartıntinn cear musat nerform successfully for the 
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system to operate. 

3. Heat exchanger must perform successfully for the 
system to operate, 

4, Salt water pump must nerform successfully for the 
heat exchanger to operate. 

5. Fresh water pump must perform successfully for the 
heat exchanger to operate, 

6. Gas turbine must perform successfully when more than 
cruising power is required. 

7. Lube oil has to be available to the diesel at all 
times. 

8, А defined shock level with a probability of pec 
occurring will cause failure of the reduction gear 
anc heat exchanger. 

Interaction effects: No interaction among the component 

failure distributions will be assumed. This is not 

necessarily the case in the real system but at the 
present state of the art, interaction effects cannot be 
handled except by using "guess" factors. Ina later 
part of the thesis, these interaction effects will be 
investigated. | 

Maintenance policies: For the propulsion plant of this 

example, no true maintenance actions will be considered 

performed. It will only be assumed that minor adjust- | 
ments for power, temperature, and fuel regulation will 
be made. Of course in most systems, maintenance actions 


will improve the reliability and this effect will be 
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covered in another section of the thesis. 
Concept of Success and Failure 

In order to determine the failure distribution of a 
component, it is зүйг Со define what is considered a 
failure and what a success in the operation of the component. 
It is not necessarily true that the onus failure of a com- 
ponent is when it T— functioning because in many 
situations, the system will also fail if a component oper- 
— certain tolerances. For instance, if a lube 
oil pump doesn't supply the oil at the required pressure 
even though the pump may still be running, the system will 
fail. It is therefore evident that the investigator must be 
careful in choosing the ground rules for success or failure 
of any test and relate them to the various modes of operation. 

To find the failure distributions relating to the com- 
ponents, a four step procedure must be followed; | 

1. Define success and failure of the entire system 
being analyzed for a particular mode of operation. 

2, Relate the definition deterinined from 1 to each 
individual components operation so that success-failure 
criteria may be ascertained for each component. 

3. Establish the type of statistical distribution which 
describes the failure phenomena. 

4. Estimate the parameters which completely define 
these distributions. 

After this is accomplished, the investigator is ready 
to begin the mathematical aspect of the analysis since all 
input data are now available. For the example of this thesis, 


t will be assumed that the above steps have been carried out 
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and the failure distributions of the components in question 
are known. The case where the out of tolerance distribution 
is independent of and different than the "catastrophic" 
failure distribution will be taken into account in one case. 
Reliability Schematic Block Diagram 

The reliability schematic block diagram sometimes called 
the reliability structure model is a picture form of the 
functional relationship of the system and the components 
under investigation. It is not a Sapam: facet of an oper- 
ational analysis, however in complex situations the model 
will provide an insight into how to best attack the problem 
from the mathematical standpoint. From this diagram it is 
also easier to recognize the weaknesses or potential weak- 
nesses of the system from a reliability sense. The relia- 
bility schematic block diagram for the thesis example is 


shown in figure (V.. 
| FIGURE V 


-= = eee š amma — s = 
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III Mathematical Procedure 
Three Basic Rules of Reliability Mathematics 

Before the mathematical analysis of the system is begun, 
a statement of the basic reliability rules of probability 
calculus will be discussed. 

The first rule states that if there are N mutually 
independent components, each having a reliability of R, (t) 
and that for satisfactory operation of the system all the 
components must function properly; then the overall system 
reliability R (t) is 


вув 


N 
в (+) = ТТ R(t) 


сув. 1 


This is known as the Product Rule of Reliability and a 
system of this type is called a serial system. 

The second rule states that if there are N mutually 
independent components each having a reliability of R, (t) and 
unreliability of Q, (t) so that 


R, (t) + Q (t) =1 
and that for satisfactory operation of the system, only one 


of these components must function properly, then the overall 


system reliability R (t) can be determined since 


sys. 
N N 
Q (6) = TT q,(t) = TT (1-2, (t)) 
сув. 1 | i-1 a 
R (t) 21 - Q (t) 


sys: sys.. 
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: N 
R (t) =1 - T T (1-R,(t)) 
sys: 1-21 Ж 


This is known as the Product Rule of Unreliabilities and 
a system of this type is called a parallel system. This is 
usually described by the term redundancy, relating to the fact 
that there are alternate components to help the system 
operate successfully in case of failure of one or more of 
the other components. The normal method for increasing 
system reliability is to add redundancy to the original system. 

The third and most important rule is based on Bayes' 
probability theorem. It states that the probability of 
system failure is the probability of system failure given 
that the component is bad times the probability the component 
is bad plus the probability of system failure given that the 
component is good times the probability the component is good. 

Q(t) = 6. (+) в. (5) x Ry(t) + O,(t)1,(t) x Q(t) 

This is called the Conditional Rule of Reliability and 

can be applied to any system including serial and parallel. 
Mathematical Analysis of the Complex System 

With these rules the analysis of the complex system can 
be begun. They will not account for every situation en- 
countered in the investigation but they will be the foundation 
for all the computational work done. 

For the purpose of simplicity, the components will be 
referred to by the number in the upper left hand side of the 
block in Figure V , i.e. 1 = diesel, 2 = gas turbine, etc. 

It is seen that blocks A, В, С, D form a serial system such 
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that if any one of these fail the system fails. If the 
reliability of each of these four blocks could be found then 


E = R,(t) -Ry(t)-R,(t) -Rp(t) I 


This will be the approach to determ ining the overall system 
reliability. 
Reduction Gear 

The first component to be examined will be the reduction 
gear. It is known that the reliability of the gear is Ry (€) 
under normal operating conditions. However, 1t 18 also known 
that some predefined shock amplitude will cause the gear to 
fail. 

Define Р as the probability that a shock of the 
defined amplitude will occur 


9. =1-Р. 


Reliability of D = Reliability of 4 given no shock 

occurs times the probability of no 
shock plus the reliability of 4 
given the shock occurs times the 
probability that shock occurs. 

Rp(t) = (Ry (4) |Р, = 0) Q, + (Ry (t)| P,P. 

(Ry (8) | 2, = 0) = Ry(t) Ry (t)| P, - 0 

therefore 


R (t) = Ry(t) Q, ET 


Heat Exchanger 
In finding the total reliability of block С, we again 


use the Bayes Theory approach of Rule 3. 





240" 
Rolt) Р. = 0\= R(t) (if rs go0d)R, (6)+R3(6) (11 7 is bad)Q, (t) 
= R3(t) (it АЕ good)R, (t) 
as В (6) (12 Т 15 bad) = O 


which leaves R.(t)(if 7 is good) to be found 
R,(t) (if 7 is good) = Ra(t) (if 8 is good)Re(t)+ 
Ra(t) (if 8 is bad )Qg(t) 


Ra(t)(it 7 is good) = На (5) (17 В 15 good }Re(t) 


as На (Е) (12 8 3s bad) = O 
but | 
Ка (5) (ат 7 and 8 are good) = Rolt B. = 0) 
^. Ro(£Y P, = 0) - RS(t)R, (t)Rg(t) 
and it is seen that the dependence of the heat exchanger on 
the fresh water and salt water pumps is the same as a 
subsystem composed of these three components in a serial 
arrangement. However, the shock hasn't been taken into 
account therefore 
Ro(t) = R3(t)R, (t)Rg(t)-Q, Tu 
Gas Turbine 
In the analysis of the blockcontaining the gas turbine 
the first unusual case of mathematical reliability prediction 
appears where the three basio rules are not sufficient 
enough to give an answer, In this case, probability theory 
must be used to formulate a new type of reliability equation, 
This is due to the fact that the gas turbine is only used 


periodically and only needs to operate properly when calied 


into use, 





са 


The reliability of the gas turbine can therefore be 
found by the following procedure; 


Define 
f(t) = failure density function of the gas turbine 
РЬ(%) = probability that the gas turbine is not іп 


use at time t 
u„(t) = distribution function of periods of non-use of 
the gas turbine 


Q(t) unreliability of the gas turbine. 


The situation can be described by the statement; a 
failure of the gas turbine at time T < t causes system 
failure by t under one of two conditions: 

1. The gas turbine is in use at time т when it fails. 

2, The gasturbine is not in use at т but is called 

into use after т and at or before t. 


Condition 1 is characterized by 


a(t) = ог (т) 1 (т)ат 


Condition 2 is 


qn (t) 


Г е (т) 0 (е-т) г (т)ат 


and 


QA (t) + 8206) Ту 


a(t) 


Since the periods of non use are considered to be 
randomly spaced, 
Defining 


N, = mean period of gas turbine non-use 


2 
the non-use density function is 
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- t/N 
=: 2 
u,(t) ти“ 
and J / É 
| - Т/М - Т/М 
U(t) = ire = e = ar . е >| 
O 2 O 
U.(t) = 1 - е ИДЕ 


Diesel Engine 
In deriving the reliability equation for block А, 


encompassing the diesel engine, another situation arises when 
the three basic reliability equations are not adequate. This 
situation is caused by the lube oil-standby lube oil pump 
combination and this combination is referred to as a standby 
system, To find the reliability of such a system 
Define 

г. (6) = failure density function of the lube oil pump 

25 (+) = failure density function of the emergency lube 

oil pump 

Since the standby lube oil pump only operates at some time 
Ey > 
functions for the system are 


when the lube oil pump fails, the true failure density 


fo (t) lube oil pump 
f (t-t) = standby lube oil pump 


As the only time the system fails is when both pumps fail, 
the system unreliability is found by Rule 2. 


A = 9(6)96(%) = г t5 (&,)at, f" fg (1-t. )àr 


and the failure distribution function of the emergency lube 


oil pump 
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Fo(t-t,) = [ £g(T-t, Jar 


R(t) = 1 - Q(t) 
consequently 


R (t) 
standby 


mean Fg(t-t,)£,(t,)aty. V 
O 


Using the same type analysis as for the heat exchanger 


Ra (t) = D VI 


Overall System Reliability 
To find the overall system reliability, equation I is 


applied 

R(t) = Ra(t) -Ry(t) -Ro(t) «my (6) г(а) 
ипеге 

849) = ыы 

Rg(t) = [1 - Qo (t)] 

Ro(t) = g(t) +R, (t) Re (t) -@, 

Ry(t) = Ry (t) ag 


Equation I(a) is applicable for the overall system 
reliability for any type failure distributions fbr the com- 
ponents t.e. exponential, normal, Weibull, etc. In most 
instances, however, the solution of such complex equations 
necessitates the use of a digital computer. For the purpose 
of easy understanding and readily obtainable results, the 
failure distributions of the components will be assumed expo- 
nential. This is the standard method in the majority of 
reliability analysis and for most complex systems not too 


unrealistic, at least for a first approximation. Тһе 
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assumption of the exponential distribution requires two basic 
Pee cuilates: | 

1. The components are not operated long enough so that 
they begin to wear out due to aging. From a reliability 
standpoint this is a logical way to operate the components 
as their reliability greatly decreases as they enter the 
wear out regime. 

2. Any failure of a component is completely random, 
i.e. it is not time dependent. 

Using the exponential distribution, the following 


equations are arrived at; 


207101 C -Ct -C-t 
Rp(t) = e ЕР сос. (е ол 


See Appendix [B ] for proof. 
= = 
-Cat C C 
R(t) =e 2 [1+ P(t) +e E Pt) EN] 
| A. xp = С 


See Appendix [B] for proof. 


-(С.+С,+Са) 5 
Ж Su ] 


Ro (t) 


-Cy t 
R(t) Q [e 


thus 
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-(с +с.+с, +С,+Са)6 
ээ ас A NE 


-C.t C C -C-t 
5 6 
sys- E | (1% ea 


ct 


times TS e 2 | 
N^ Co: N, ^ C5 
Quantitative Analysis 
In order to more graphically show the value of the relia- 
bility analysis, representative values for the component 
instant failure rates (po (t) = С.) will be chosen. The 
overall system reliability де а Amen of time can then be 


evaluated, 


The values chosen are; 


Су = .0002 np”? 
-1 
Cp = .0003 hr 
Сз = .00006 np”? 
Cy = .00004 hp”? 
21 
СБ = ,000À hr 
Ed 
Ср = .0005 hr 
21 
С, - .0005 hr 
= 
Cg = 20005 hr 
о. = .98 
ts) = „ВО 
No = 200 hr 


Also, assume the wear-out of only one component is taken 


into account. This component's failure due to wear-out shall 





-36- 


be described by a normal distribution with a mean life time 


of 1200 hrs and a varience of 400 hrs. 

Figure VI shows R (t) vs t for the cases with and with- 
out one component кте. дует та ий 

For the ригрове of understanding the individual com- 
ponents pélationshin to the system as a whole, each of four 
component failure distributions have been varied to find its 


total effect on the system. These four components were chosen 


due to the different nature of their operation; 


Component Time Operation 
Diesel continuous, dependent upon L. O. system 


Gas turbine periodic 
L. O. pump continuous with standby system 


Е. W. pump continuous, nondependent 
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IV Preventive Maintenance: 

In most of the reliability work done in the past, little 
interest has been given to the effect of preventive mainte- 
nance . The reason for this is twofold due to the fact that 
most of the reliability study has been applied to the missile 
and electronic fields; 

1. Іп the missile field the operation is usually "one 
shot" i.e. the component or system having once begun its 
operation cannot be maintained, and after terminating this 
operation it is unusable. 

2, TIn the electronic field, the failure distributions 
of the components are considered exponential thus making 
preventive maintenance useless. This is so because the 
instant failure rate is independent of time and therefore 
replacing the component at some specific time when it is 
still operating successfully has no advantageous effect from 
a reliability standpoint. 

These two conditions do not necessarily hold however, 
for a marine propulsion plant. Such plants definitely are 
not "one-mission" devices and it is doubtful if the exponen- 
tial distribution will be applicable to the mechanical com- 
ponents found in such systems. It is proven by R. F. Drenick 
[4] that if the instant failure rate is constant or decreases 
with time, then preventive maintenance is not beneficial. Іп 
mechanical components, it. has been determined that the instant 
failure rate normally increases with time. In this case pre- 


ventive maintenance can be worthwhile and should be considered 
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in any operational analysis of a marine propulsion system. 

A maintenance investigation should provide two basic 
parameters; how valuable to system performance iS preventive 
maintenance, and at what time intervals should it be per- 
formed, For the purpose of analysis two different replacement 
policies will be explored. The mathematical analysis is from 


the work of R. Barlow and L. Hunter [5]. 


Бой су І: 

After every ES hours of continuing operation without 
failure, preventive maintenance is performed. If there is a 
failure before to» the maintenance is performed at this 


time, 21 and the next maintenance is scheduled at time 


t1 


+ ®. 
Defining 
No = expected fractional amount of time the system 
is operating as t — ° 
the two parameters that ed the condition of the problem 
are п. апа +. The criterion of optimality will be the 
maximization of T, ° 
Те = expected time needed to perform repair to the 
system after a failure 


= 
II 


expected time needed to perform a scheduled 
preventive maintenance action. 

By making the assumption that after each repair action, the 
system is as good as new, the following equations are appli- 


cable for optimal results; 
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5 T 
p(t,) J ° R(t)at - F(t) = m- 


O e S 


1 
T ES IFT -T Jolt F 
1 со 1+ T. T, p A 


if To < To no preventive maintenance should be made and 


Tp 


Tho 7 TRUE 


Policy II: 

In this case, preventive maintenance is performed on the 
system after it has been operating for a period of t? hours 
regardless of the number of intervening failures. If there 
is a failure before t° hours have elapsed only minimal 
repair is made which in effect puts the system back in oper- 
ating condition but does not change the basic system failure 
rate. 

Using the same definitions for m, and T, and defining; 

В = expected time needed to perform minimal repair 


to the system 
the following equations are determined 


Ls 23 
f tp'(t)àt = т 
O m 

_ 1 
Nov d 


O 
14T, p(t ) 
The question naturally arises as to which policy should 
be used for a particular situation. This question can be 


resolved by equating 14 = Mo for a given Т, which is 
known and the same for both policies. By substitution of 
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values for Те, Іп can be determined. By graphing 
Те vs In for a given Ta» which policy to use becomes 
evident. Figure XI shows such a graph for a Weibull 
distribution. 
A Maintenance Policy by Dynamic Programming 

Assume that a preventive maintenance policy is to be 
determined for a system having an instant failure rate that 
is a linear function of its operating time. The system is 
to be operated in discrete time intervals of length бо. 
ТР the system is repaired between operating periods the cost 
ses ateq with the repair is С * If the system fails 
during an operating interval it cost Ce for the repair and 


the system cannot be operated again till the next interval. 


Define 
С = cost of seheduled maintenance 
Cp = cost of "in service" repair 
P, = Probability of failure in interval m+ 1 


given it has not failed in m intervals 
since the last maintenance 
ғ (1) = least expected cost of making n more 
| intervals given m intervals have been 
made. 

There are two policies that can be followed, either 
repair the system at m or do not repair it. By the 
methods of dynamic programming covered by R. Bellman [6 ] 
the following dynamic programming equation is developed; 





IT 


eine 


R: CtP [ct _,(0)] + (1-Р,)1, 104); repair 
ғ (т) - тіп 


NR: PaL Oe fa 1 (1 + (1-Р 27, Q4 (001); non-repair 


Се > С for logical results otherwise preventive mainte- 


É 
nance is not worthwhile. 

To show mathematically the effect of the failure repair 
cost, Ce: and the scheduled maintenance cost, С.» on the 
maintenance policy to be followed, the following problem will 
be postulated. 


A system has a failure distribution curve that can be 


described by the Weibull approximation with the following 


parameters ; 
y = location parameter = O 
B = shape parameter = 2 
a = scale parameter = 57 x 10" hrs 
о2/8 = characteristic lifetime = T = 755 hr 
M(t) = N(t) = interval time = 100 hr 
therefore 
alt 
сп 
+2 
= —— 


R(t) = e 210048 ) 
The question to be answered is how does the maintenance 


policy vary as a function of m, n, Са» апа С г Figures 


XII through XVI give the results of this investigation. 
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V Interaction of Component Failure Distributions 

In the previous sections of the thesis, the overall 
reliability of a hypothetical complex marine propulsion 
plant has been derived by the proper mathematical combi- 
nation of individual component reliabilities. In actual 
situations however, it has been found that although this 
method is mathematically correct, it does not yield the 
actual reliability observed in practice. Intuitively, it 
might at first glance appear that this poor correlation is 
because the model is not a good functional representation 
of the real system. Although this is a possible reason, 
further analysis may show that mathematically all work is 
accurate. Where, then, is the mistake in the analysis? 
One of the major possible explanations and one that although 
recognized is not included in present reliability studies, 
is that when components are integrated into the systen, 
there is an interaction among them which can change their 
failure distribution. Because of this, the failure distri- 
bution observed in the single component life test is not 
applicable and must be modified. The problem of finding what 
type of function this modification factor is, will in actual 
practice be very difficult. However, if reliability investi- 
gations are to be accurate, component failure rate inter- 
action has to be taken into account as it may have a 
pronounced effect on system reliability, preventive mainte- 
nance policies, and system availability calculations. It is 
therefore imperative that laboratory life tests of components 


suspected of interaction and data assimilation of real system 
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operation be analyzed so that interaction functions can be 
determined. | 

For the purpose of illustration, two possible situations 
where interaction may play a substantial role in an investi- 
gation are given. | 

a, The failure distribution of a diesel cylinder liner 

will be affected by the wear of the biston rings, 
the wear of the connecting rod bearings, and the 
fuel quality. 

b. The failure distribution of a reduction gear will 

be dependent upon the wear of the shaft bearings and 
coupling operation. 

Since no data is available on interaction, for the 
purpose of investigation only hypothetical interaction 
functions can be postulated. In the following sections such 
assumed functions are examined. 

Effect of Interaction on Component Reliability 

The time dependent reliability of a component is de- 

scribed by 
R(t) = expl- Г p(t) dt] 


which is a function of a single variable, p(t) = the instant 
failure rate. An equation illustrating the шалт factors that 
influence component failure rate is characterized by 

p(t) = f(a, B, 6, t) 


derating factor 


= environmental stress factor 


interaction factor 


+t © D R 
N 


operating time 





u. 


where 

a= accounts for the fact that the component can operate 
at different outputs 

B = is a function of the operational environment the 
component experiences, i.e. temperature, vibration, 
and acceleration stress levels 

b = accounts for the effect of interaction among the 
component failure distributions 


t = time the component is in operation 


It is normally possible to find the instant failure rate 
as a function of derating factor, stress factor, and time 
from experimental life testing of the component. An example 


of such a function can be shown in a diagram for a component 


that follows the exponential failure distribution. 
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Because reliabllity engineering is relatively new, little 
experimental data is available to find such functional re- 
lationships and what is available only applies to the expo- 
пепбїа1 failure distribution. It is valid to assume that the 
same type of relationships or equations can be found for other 
types of components and failure dsitributions. The problem 
still presents itself 57 how can the interaction factor Бе 
handled? Since no work has been done on this matter, only 
conjectural relationships can be assumed for the purposes of 


investigation. Such relationships will therefore be assumed 


and their properties analyzed. 


Define 
A(t) = component instant failure rate without taking 
| into account this interaction factor, %. 
therefore | | 
p(t) = Ғ(>(6), %) 


It is this function that will be hypothesized. In all 
cases A(t) will be considered a constant for ease of mathe- 
matical computation. The effect of the interaction factor 
which is determined in this exponential case will apply in 
effect to other failure distributions. 


Investigation 1: = additive and constant 


p(t) = À + Z C 


и 1 


R(t) = exp[-(A + 2 C, )t] 
| RS 


Те = 1/(3 + 4 C4) 








таз 


Investigation 2: Ó = additive and time dependent 


op(t) = À + 5 C.t 


R(t) » exp[-^t] expl- >> ct" 
| n 


= әр {Е EE 0. D 


i 
See Appendix[Blfor proof. П 


Investigation 3: Ó = multiplicative and time dependent 


p(t) = (ct) À 
R(t) = exp[- ЗА t°] 
Т 


- МЕН, 
| (cA) 
Investigation 4: = additive and a function of a constant 


and a time dependent term 


e(t) БА + С. + Cot 


from correlation with investigation 2 


R(t) = exp[- (C4 )t] expl- > 664) 
2 
E A+C (12 егер С | 
Te exp[ | ] | BC ( er [ >? 


It can be seen that these values for R(t) and Те 
сап differ greatly from the values obtained for the standard 
exponential failure distribution, i.e. 

R(t) = exp[-at] 

E X 
Figure XVII shows the relative importance of > С. for 


n 
investigation 2 where A= .Ol. 
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FIGURE XIII 


Preventive Maintenance Policy by 
Dynamic Programming for Ce = 10 
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FIGURE XV a 


Preventive Maintenance Policy by 
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FIGURE XVI 


Preventive Maintenance Policy by 
Dynamic Programming for ratio 
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DISCUSSION OF RESULTS 


Reliability Analysis 
In Figure VI the overall system reliability as a function 


of time is shown. The plot graphically displays what degree 
of reliability can be expected from the system assuming; 
а. the values for the instant failure rates specified 
previously. 
b. the components are good as пеи ат ї = 0. 
с. no preventive maintenance is performed on the system 
during the operating cycle. 
It can be determined that this overall reliability curve 
for the system consisting of different components with unlike 
operating conditions can be approximated by; 


-Ct 


R (t) =e 


sys- 


where C = .00198 
Other results that can be obtained from Figure VI are; 
l. If the instant failure rates of all the components 
were decreased by 50%, the overall system reliability will be 


increased greatly as t increases, i.e. 


t = 200; Куз is increased by 20.2% 
t = 500; Вуз is increased by 61.2% 
t = 800; Roys is increased by 116.8% 


2. If one component has wear-out characteristics (fresh 
water pump) such that it should be included in the reliability 
analysis, the overall system reliability is adversely affected, 


This effect becomes especially prominent as the system 
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operating time approaches the mean life time of the component 
and is also dependent upon the varience of the mean life time. 


Figure VI shows that at t = 800, R із decreased by 15.7% 


sys 
due to wear-out. The equation that is applicable when investi- 
gating a system or component which exhibits random and wear- 
out failures during the operating cycle is; 


R, (t) = Ва) Ry (t) 


defining R (t) = overall reliability of component i 


R, (t) = reliability of component i exhibiting 
и 
only wear-out failures 
R, (t) = reliability of component i exhibiting 
qe 


only random failures 

Since the wear-out reliability factor is multiplicative, 
the system reliability will decrease greatly with time as the 
number of components experiencing wear-out increases. 

3. Figures VII through X illustrate the relative 
importance of each component to overall system operation 
from a reliability viewpoint, These figures show percentage- 
wise how the decrease in a specific components instant failure 


rate will affect the system by using the formula; 


№ s (with decreased p(t))-R, „(04 (5)-с.) 
2 Bsys P4 9794. 

It is observed that a decrease in the fresh water pumps 
failure rate by 50% increases the overall system reliability 
by 22% for an operating cycle of 800 hours while the same 


decrease for the diesel engines failure rate only increases 








дави 


the reliability by 7.8%. Therefore if the system reliability 
needs to be improved, strictly from a reliability standpoint 
of the four components investigated the fresh water pump per- 
formance should be improved first and then the lube oil pump, 
gas turbine, and diesel respectively. An investigation like 
this immediately indicates where in the system yimprovements 
should be made to increase overall reliability. Conversely, 
it also reveals what increase can be expected in system relia- 
bility for an improvement in a specific components instant 
failure rate. Consequently, graphs of the type shown in 
Figures VI through X, which can be determined for any system 
if the instant failure rates of its components are known, 
provide an intimate insight into system behavior and provide 
valuable information for decisions regarding engineering 
improvements. 
Preventive Maintenance 

Figure XI shows the relative merit of using either of 
the two policies discussed on pages 38 -409 of this 
thesis. This graph is representative of the preventive 
maintenance situation for a system having a failure rate 
that increases with time. As an example of how to use the 
figure, assuming T. = 7 hours and the expected value of 
Те = 1 hours then policy I is applicable unless Басе 
hours, Similarly, for any value of To нава is а согге- 
sponding value of In which determines the boundary for 
choosing the optimum policy to follow. It will generally 
be found. that policy II is most applicable to complex systems 
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which have many subsystems liable to fail and where keeping 
records of times to failure for each subsystem becomes 
cumbersome, It is interesting to note that if costs of 
failure are substituted for the times of repair, then these 
same equations minimize the expected cost of operating the 
system. 

The method of determining preventive maintenance 
scheduling by dynamic programming is more suitable to a 
marine propulsion plant than the previously discussed method. 
The reason for this is that with operating schedules to main- 
tain it is probable that optimum preventive maintenance times 
found will not coincide with the ships schedule. This method 
also provides a continuous picture of what policy to take re- 
garding preventive maintenance over all combinations of time 
intervals (trips) completed and time intervals still to be 
completed. Б рев XII through XVI illustrate the results of 
such an analysis made for a system having the Weibull failure 
distribution described in the Procedure section. From these 
graphs the expected cost of maintenance can be minimized for 
any operational situation of this system if the cost of 
failure and of maintenance are known. These figures are only 
applicable for m, n> O as if m= O, the system has just 
been repaired and it does no good to repair it again, and if 
n = 0, there are no more trips to make so the expected cost 
= 5 

Although the graphs are self-explanatory, one interesting 


observation can be made which is not intuitively obvious nor 
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would normally be expected if the investigation had not been 

made, Take as an example, Figure XIII for Ce = 10 and the 

curve representing С, - 2. It is seen that if 8 trips have 

been made (m = 8) and 2 more are planned (n = 2), the policy 

should be to пал the system. However, with the same past 

history and if 3 more trips are necessary (n = 3), the system 
should not be repaired. The explanation of ме” 

l. at large m the probability of failure in n inter- 
vals is relatively high and for small n, repair is the best 
eier since in effect it starts the system off at m= 0, 

2, as n increases, this effect of starting at m= 0 
becomes less important and the cost of preventive maintenance 
is the dominant factor. The reason is that at large n, even 
if maintenance is performed, there is still a significant 
probability of failure in the following 1-1 intervals which 
causes a greater expected cost than if no repair is made and 
a high probability of system failure is accepted in the next 
n intervals. 

Analysis of Figures XII through XV provides the results 
for the determination of Figure XVI. Figure XVI shows that 
policy determination is solely dependent upon failure distri- 


bution, length of interval, and the ratio of C to С, 


T D 

This fact lessens greatly the number of calculations that 
have to be made in a system analysis since assuming the 
system failure distribution and interval time are known, only 
one graph has to be derived instead of one for each possible 


value of Cre As the values of Ce and C, are bound to 
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vary with time due to the normal trend of increasing repair 
and material costs, the advantage of the single graph repre- 
senting the ratio is evident. 
Interacting Components 

Although little if any investigation has been made into 
the problem of the interaction component failure distri- 
butions, Figure XVII illustrates that in certain cases this 
effect can be very important. For instance, it is seen that 
if p(t) = .01 + 4(102)t the mean time to failure is 
а 20%. On the e hand for t =“ sf, this inter- 
action effect only decreases reliability by less than 1%. 
Consequently although the effects of interaction may not 
affeet the basic calculations for reliability significantly 
and therefore may be difficult to observe and determine, they 
definitely can become ге in the overall system 
analysis through Te which is used in many operations analy- 
sis calculations not discussed in this thesis, i.e. renewal 
theory. It therefore becomes imperative that investigations 
should be made determining any interaction present in a 


system if worthwhile results are to be obtained from an 


analysis of the system. 
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CONCLUSIONS 


Tn addition to the specific conclusions stated in the 
Discussion of Results, the following general conclusions can 
be drawn. 

1. The operational analysis of a marine propulsion 
plant will provide results not otherwise available which 
will be beneficial in making decisions regarding; 

а. preliminary design 

b. operational policies 

с. maintenance actions 

а, system improvement and redesign 

2. In such an analysis the definition and formulation 
of the objectives of the investigation determine all suc- 
ceeding stages of the analysis and the utility of the results, 

3. Statistical data are a necessary prerequisite to any 
operational analysis. The definition of the problem governs 
the areas of investigation for the statistical tests. 

4. Reliability as a function of time can be character- 
ized by the single parameter p(t), instant failure rate. 
The effect of an increase or ее in this parameter is 
amplified as time increases. 

| 2% The wear-out failure distributions of the individu- 
al system components should be included in the analysis of 

a mechanical system as they can greatly affect the overall 
| system reliability. This consequence becomes especially 
prominent as t approaches the MTBF of a component, 


6. From an operational analysis of a system, the 
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pelationship of each component's performance to overall 
system operation can be determined. Such information will 
be of considerable significance in the determination of 
policies regarding system design and improvement. 

7. Preventive maintenance is only advantageous for 
devices having an instant failure rate that increases with 
time. 

8. For this type device there is an optimal time for 
the performance of preventive maintenance which will maxi- 
mize availability. An optimal time may also be found which 
will minimize expected cost of operation. 

9. Dynamic programming can be utilized in determining 
an optimum maintenance policy for systems operating in 
discrete time intervals. This policy will give a complete 
picture at any time in the systems operating cycle of what 
maintenance action to follow in order to minimize the 
expected cost of repair. 

lO. The effect of interacting component failure distri- 
butions normally neglected in operational analyses and 
difficult to determine quantitatively can considerably 
affect the results of the investigation and should therefore 
be included. 
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RECOMMENDATIONS 


1. The concept of the operational analysis of marine 
propulsion plants should be accepted by naval engineers as 
an important tool in the evaluation of existing and proposed 
marine systems. 

2. A program of component testing and data assimilation 
should be initiated immediately for the purpose of determi- 
ning representative failure distributions of mechanical 
devices. 

3. Accepting the fact that corroborating statistical 
data are minimal, systems analyses should be carried out on 
existing systems for the object of developing and refining 
techniques and gaining facility in the procedures of such an 
investigation. 

4. Further examination should be made into the question 
of how out of tolerence i.e. wear-out failures are to be 
defined and handled. 

5. A detailed investigation should be made into the 
problem of failure distribution dependence. This should be 
explored both from a theoretically mathematical and an 
experimental viewpoint. 

6. The applicability of the method of dynamic program- 
ming to system analyses questions should be investigated. 
Possible applications other than for maintenance policies 
might Те; 

a. optimal reliability apportionment among the 


components, 
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optimal arrangement of components under con- 
straints for maximum system reliability, 
optimal choice of components to minimize the 
deleterious effect of failure distribution 


interaction. 
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APPENDIX A 
Definitions 
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Failure: а detected cessation of ability to perform a speci- 
fied function or functions within previously established limits. 
Independent Failure: those Amen failures which occur or 
can occur without being related to the malfunctioning of associ- 
ated components. 

Redundancy: the existence of more than one means for ac- 
complishing a given task where all means must fail before 
there is an over-all failure to the system. 

Availability: the fraction the total desired operating time 
that a system is actually operable. 

Mean Time Before Failure (MTBF): the mean or average time 
between successive failures of a component 

Cumulative Probability of Failure (Q(t), F(t)): probability 
that the life of a component is less than t. 

Probability Density Function of Failure (f(t)): probability 
that the life of a component will be between + авт t + dt. 
Instant Failure Rate (p(t)): probability that a component 
will fail during the time interval t to t + dt conditional 
upon it surviving up to time t, 

Reliability R(t): probability that a component will operate 
within specified limits for the time, t, and operating con- 
ditions imposed. 
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APPENDIX B 
Details of Procedure 
I Periodically used component 
ІІ Standby system 
III Interaction 


TV Computer program 
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I, Determination of R(t) for a component used periodically 


assuming an exponential failure distribution. 


Р(+) = /® [1-Р(т)]2(т)ат + /® Р(т)[6(%-т)]г(т)ат 
O O 








P(t) = P 
x x (t=) 
0(%-т) -1-е N 
f(t) - С ЕБ" 
Q (T) = ie [1-P]C е Стат = [1581 [1-2 #1 
| O 
21720 P3. «СЪ + P «СЪ 
1 
= att 
Q(t) = ie Р[ 1-е х o o CT aq 
i O 
t 
= P 1 C + p(_£ N 
е "те Tw 
„+ 
F(t) 21 - e"C*[3 4- P(7 С )] + P(+ Ç )ë N 
| N к-с 
R(t) = 1 - F(t) 
t 
N 


R(t) = + A) - r(e 
| NC С. ў - С. 


II. Determination of R(t) for a standby system assuming 
exponential failure distributions. 
p(t) = С. for main component 
p(t) = C, for standby component 
_ pt 
For a standby system Q(t) = ) Fo(t-t,)1, (t, Jdt, 
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ОЕ) 
Bic = Ғы = Пе “° la 
253211 1 2 
t t 
1 1 
"a -C, (t-t,) 
t t -С2(%-%)) | -С|% 
Q(t) = J Fo(t-t,)1, (t)dt, = J [l-e ]C,e dt, 
-C.t -C.t (6.-C. )t 
Б E cie * lat, - Qe 2 [ina e * 1. lat 
1 1 1 
О O 
t 
-C. t С -0-5 . 08--0-)5 
1 1 2 DO VI 
=1l-e - ele “ )[e =] 
281 O 
-et С Zeit -С85 
un te 52 oe 1 -е < ] 
R(t) = 1 - Q(t) 
С -с 
1 1 2 
R(t) = e [1 + ] - e 
С2-С, 


ніні» 


Determination of R(t) and T, for p(t) = A +C,t 


f 1 


li 


p(t) = à + СЪ 


R(t) = expl -fË (a + с, т)ат! 
O 


= expl -ft лат - K Сут ат] 
O O 


= exp[-Aat - С. i т ат! 
O 


С 
= exp[-At - 5 65] 


2 


1 
- 5 C,t 
de LEM mE 


T. = f° R(t)dt 
O 








ST = 


со 1 2 
= J exp[ -(^t % 2 6,% )] at 


From C.R.C. Standard Mathematical Tables [7] 
oo 2296 | _ 4887 
J ехр[-а“х“] ах = 5= Ут = — 


it can be seen that 


2 2 
d 2 À À 
C,t + № = 5 с, Her 5+ 2. E 
moy С, C, C, | 
2 2 
1 À À 
= C.[(t +—) - (—) ] 
2 1 С.“ C," 
2 
= 5 C(t + оны кс, 


therefore 


R(t) = explo expl- 5 C (t + с) 3 
therefore 


2 
= f” R(t)dt = ача Г expl- $ C, (t + 2-) 146 
O | 3. 


let x = t ++ 6-0 x = 
i 


a» 


8 


dx = dt t = x = оо 


со 2 со 
J exp[- 5 C, (t +7) 1456 = ! exp [ - 5 с, (к ах 


Q 


which equals (dropping the subscript) 
A 
Dm ехр( x2]ax - f? exp[Z x°] ax 
O O 


from C.R.C. tables 


со со T 
J exp[5 x ]àx N 
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now let $ = С х x = 0 t = 0 
at = VS ах x = à ue X 
JOE 
therefore 


A 


A 
re ехр! хе ах = 29 (D) E мес exp[ -t^]àt 


O O 


and À 


2 pe 2 À 
Мт о E | “ Л20 


thus 
P exp[2 х“)ах = JE [1 - ert (——)] 
and 


2 = 
Th = expl sg] IT. (1 - ere) 
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IV Main Body of Computer Program for Solution of Equation 


100 


104 


102 


103 


104 


105 


106 


107 
108 


109 


(FORTRAN) 

pò 114 SH =1, In 
JJ = O 

DO 114 M = 1, М1 


IF(N-1) 100, 100, 101 
Al = 0, 


GO TO 102 
2-М-1 

Al - cosmi(J,1) 
A2 = CASREP + Al 
A3 = RATE(M) * A2 
IF (N-1) 103, 103, 104 
Ak = 0. 

GO TO 105 
К-М-1 

A4 = COST1(J,K) 
A5 = 1.-RATE(M) 
A6 = A5 * AL | 


COST1(N,M) = АЗ + Аб 
TF(N-1) 106, 106, 107 
Bi = O. 


GO TO 108 
Bl = COST2(J,1) 


B2 = SKDREP + RATE(1) * (CASREP + B1) 
IF(N-1) 109, 109, 110 
B3 = 0. 








110 
жы 


да 


113 
114 


- ВЕ 


со то 111 

B3 = COST2(J,2) 
B4 = 1.-RATE(1) 
В5 = Ві * ВЗ 


COST2(N,M) = B5 + B2 

IF(JJ) 112, 112, 113 

баг! 

РЕТМТ 20 

РВЇМТ 21, М 

L=M-1 

PRINT 22, L, COST1(N,M), L, COST2(N,M) 
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